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Multiple  Scatter  Cross  Sections  for 
Anisotropic  Kolmogoroff  Turbulence 


I.  INTRODUCTION 

1* 

In  a  recent  report  on  backscattering  from  anisotropic  media,  we  developed 

a  method  for  calculating  multiple  scatter  cross  sections  based  on  DeWolf's 

o 

cumulative  forward-scatter  single -backscatter  approximation.  T ne  method  con¬ 

sisted  of  obtaining  an  expression  for  the  scattered  wave  in  terms  of  the  Green's 
function  for  a  wave  propagating  in  an  arbitrary  random  medium.  From  this,  an 
expression  for  the  power  scattered  into  an  arbitrary  direction  was  obtained.  The 
scattered  power  obtained  from  averaging  over  the  random  properties  of  the  medium 
was  then  expressed  as  the  Fourier  transform  of  a  cumulant  expansion  of  the  corre¬ 
lation  functions  for  the  random  medium.  To  evaluate  the  Fourier  transforms,  the 
cumulant  series  was  further  expanded  as  a  functional  Taylor  series  in  terms  of 
the  projected  correlation  functions  (structure  functions)  of  the  medium.  Fourier 
transforms  of  products  of  arbitrary  powers  of  correlation  functions  projected 
along  the  incoming  and  scatter  directions  were  then  expressed  as  convolution 


(Received  for  publication  7  February  1985) 

♦This  report  will  hereinafter  be  referred  to  as  "I." 

1.  Yukon,  S.  P.  (1982)  Backscatter  From  Anisotropic  Random  Media, 

RADC-TR-82  -287. 

2.  DeWolf,  P.A.  (1971)  Electromagnetic  reflection  from  an  extended  turbulent 

medium.  Cumulative  forward-scatter  single -backscatter  approximation, 
IEEE  Trans.  Antennas  Propag.  AP  19:254. 
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integrals.  For  Gaussian  correlation  functions,  the  required  integrations  could  be 
evaluated  analytically.  This  calculation  was  carried  out  in  I. 


3  4  5 

Many  measurements  of  ionospheric  irregularity  fluctuation  spectra  *  *  sup¬ 


port  some  form  of  power  law  for  the  fluctuation  spectrum  and  for  such  power  law 
spectra,  the  convolution  integrals  cannot,  in  general,  be  evaluated  analytically. 
Therefore,  one  objective  of  this  report  will  be  to  evaluate  multiple  scattering 
cross  sections  for  fluctuation  spectra  that  cannot  be  convolved  analytically  and 
for  the  three-dimensional  Kolmogoroff  spectrum  in  particular.  We  will  also 
extend  the  cross  section  calculation  in  I  to  include  scattering  into  forward  direc¬ 
tions  and  scattering  from  slabs  containing  irregularities  whose  symmetry  axes 
are  oriented  at  an  arbitrary  angle  with  respect  to  the  slab  face. 


2.  CROSS  SECTION  FOR  FORWARD  SCATTERING 


The  only  formal  difference  between  the  calculation  of  scattering  into  forward 
directions  and  scattering  into  backward  direction  lies  in  changing  the  path  length 
for  the  scattered  wave  as  it  appears  in  Eq.  1-45.  Thus,  the  expression  for  the 
backward-path  length  x/?f  should  be  replaced  by  the  path  length  (L-x)/3f  which 
represents  the  path  length  from  the  scatter  site  on  the  plane  located  x  units  of 
length  from  the  front  of  the  slab  to  the  back  of  the  slab  at  L.  Integrating  over  the 
volume  of  the  slab  yields  the  analog  of  Eq.  (1-46). 


A  Jdxe 


xjff-CD.U?)-!]  +  [L-x]/3f  Df(A?)-l] 


■[AL] 


L/S.fD^A?)-!]  [Df  (Ar)-l] 


(1) 


L/3.[D.(A?)-1]  -  L/Sf[Df(A?)-l] 


3.  Erukhimov,  L.  M. ,  Kosolapenko,  V.  I. ,  Lerner,  A.M. ,  and  Myasnikov,  E.  N. . 

(1982)  Form  of  the  inhomogeneity  spectrum  of  the  high-latitude  ionosphere 
in  the  direction  of  the  geomagnetic  field.  Radio  Phys.  Quart.  Electronics 
24:350. 

4.  Kelley,  M.  C. ,  Livingston,  R.C. ,  Rino,  C.  L. ,  and  Tsunoda,  R.T.  (1982)  The 

vertical  wave  number  spectrum  of  topside  equatorial  spread  F:  Estimates  of 
backscatter  levels  and  irrp  lications  for  a  unified  theory.  J.  Geophys.  Res. 
87(A7):5217. 

5.  Ossakow,  S.  L.  (1979)  Ionospheric  irregularities.  Rev.  Geophys.  Space  Phys. 

17:521.  - - 


To  expand  this  as  a  double  functional  Taylor  series  in  D^(4r)  and  Df  (4r),  we 
first  define  Vj  =  L/?j  D.  (4r),  c.  =  L/3.  and  an  equivalent  vf.  cf  with  the  sub 
script  f  replacing  i.  The  Taylor  expansion  is  then  given  by 


vi-  ci 


vf-  cf 


[AL]  n?0  [’*  %  *  V(  ['“1  -  °i>  -  <*f  -  «t>] 


i*Vf  =  0 


(2) 


Expanding  the  bracket  using  the  binomial  theorem  yields 


[v] 


“ci  -cf 
e  -  e 

-<ci  -  V 


-rv 


•  n  1  n! 

£  E  ITT 

_ i  •  —  rs  11  ■ 


n=l  j=0 


(n- j) !  j! 


a-i  9l  /eV  c»  -  -  CA 

av^"3  avj1  y[(v.  -  Cjj-fvj-cj)]  / 


vi-vf  = 0 


(3) 


Writing  v_=  v^,  -  vf  and  c_  =  c.  -  cf,  the  derivative  term  maybe  rewritten  as 


Using  the  result  of  Eq.  (1-53),  this  may  be  re-expressed  as 


3  • 


3 

E 

k=0  (j-k)!k! 


-c.  ,  y  (n-j+k+1,  c.  -  cf) 

e  1  (-IT  1  1 


(c.  -  cf)n'3+k+1  * 


(6) 
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Comparing  Eqs.  (3)  and  (5)  with  Eq.  (1-64)  shows  that  scattering  into  the  forward 
direction  will  be  given  by  Eq.  (1-54) 

[  .  .-**•  ] 


with 


(6) 


L  (^-^f)L  J 


replacing 


1  -  e 


•(/3i  +  /3f)  L 


and 


L  (/3.  +  /9f)L 


±  j!  -0J-,  . 

y(n  +  l.j;  (/3.  -  /3f)L)  =  E  -  e  (-l)Ky(n-j+l.  (fl.-/3.)L) 

1  1  /j  l.v  m-  •  1  *■ 


(7) 


k=0  (j-k)  !k! 


j-k  (8) 


replacing 


V(n+  1.  (/3t  +  /9f)L). 


(9) 


In  calculating  cross  sections  thus  far,  we  have  assumed  that  kQLo  »  1 
where  Lq  is  the  smaller  of  the  inertial  scale  lengths  (L,,  ,  LA  ),  and  that  the 
cumulative  forward-scatter  single -backscatter  (CFSB)  approximation  holds.  This 
is  a  good  approximation  for  scattering  into  backward  directions.  However,  for 
forward  directions,  it  is  possible  that  double  and  higher  order  scattering  with 
angles  larger  than  6C  =  l/kQLo  may  contribute  appreciably  to  the  cross  section. 
DeWolf  has  calculated  that,  for  a  single  back  scattering,  the  ratio  of  scattered 
intensity  outside  a  cone  of  half  angle  9-q  to  the  intensity  inside  the  cone 
1(6  >  6^)1 1(6  <  6d)  is  approximately  given  by  (2kQLo  sin  &D/ 2)-®^.  If  this 
ratio  is  still  small  for  forward  directions,  the  single  scatter  CSFB  approximation 

should  be  adequate.  As  an  example,  if  we  require  this  ratio  to  be  no  greater  than 

-2  -1 
10  ,  this  would  imply,  for  the  choices  of  K  =  8jr.  L  =2.0  and  6  =  (K  L  )  = 

o  o  coo 

1. 14°  ,  that  =  18. 1°  .  Thus  we  should  be  reasonably  safe  in  pressing  the  CFSB 
into  within,  say,  36°  of  the  forward  scattering  direction  for  this  choice  of  para¬ 
meters. 
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Comparing  Eqs.  (3)  and  (5)  with  Eq.  (1-54)  shows  that  scattering  into  the  forward 
direction  will  be  given  by  Eq.  (1-54) 


with 


-/SfL  -0£L  • 
e  *  e _ 


(6) 


replacing 

-<A  +  0f)  L' 

1  -  e _ 

.  <0.+0f)L 


(7) 


and 


7(n  +  l.j;  03.  -  pt)  L) 


j! 


k=0  (j-k)!k! 


*/3fL  k 

e  (-1)  ytn-j+l.t^-^jL) 


J'k 


(8) 


replacing 

V(n+ 1,  (/3.  + /9f)L).  (9) 

In  calculating  cross  sections  thus  far,  we  have  assumed  that  kQLo  »  1 

where  Lq  is  the  smaller  of  the  inertial  scale  lengths  (L„  ,  Lj.  ),  and  that  the 

cumulative  forward-scatter  single -backscatter  (CFSB)  approximation  holds.  This 

is  a  good  approximation  for  scattering  into  backward  directions.  However,  for 

forward  directions,  it  is  possible  that  double  and  higher  order  scattering  with 

angles  larger  than  0C  =  l/kQLo  may  contribute  appreciably  to  the  cross  section. 

DeWolf  has  calculated  that,  for  a  single  back  scattering,  the  ratio  of  scattered 

intensity  outside  a  cone  of  half  angle  6 D  to  the  intensity  inside  the  cone 

1(0  >  0D)/I(0  <  0D)  is  approximately  given  by  (2kQLQ  sin  0D/2)~5^.  If  this 

ratio  is  still  small  for  forward  directions,  the  single  scatter  CSFB  approximation 

should  be  adequate.  As  an  example,  if  we  require  this  ratio  to  be  no  greater  than 
•2  -1 
10  ,  this  would  imply,  for  the  choices  of  K  =  87T,  L  =  2.0  and  6  =  (K  L  )  = 

O  O  0  o  o 

1. 14°  ,  that  0-p  =  18. 1°  .  Thus  we  should  be  reasonably  safe  in  pressing  the  CFSB 
into  within,  say,  36°  of  the  forward  scattering  direction  for  this  choice  of  para¬ 
meters. 
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3.  KOLMOGOROFF  SPECTRUM  OF  TURBULENCE 


The  choice  of  a  Gaussian  spectrum  of  turbulence  in  I  fnade  it  possible  for 
all  convolution  integrations,  such  as  that  in  Eq.  (1.61)  to  be  evaluated  analyti¬ 
cally.  For  power  law  spectra  afid  spectra  having  more  complex  functional  forms, 
numerical  integration  will  generally  be  required.  For  these  cases,  a  further 
distinction  may  be  made  depending  on  whether  or  not  the  projected  correlation 

functions  to  be  convolved  in  the  k  -k  plane  can  be  transformed  into  functions 

x  y 

having  circular  symmetry.  For  functions  that  can  be  represented  this  way,  the 

cross  section  calculation  may  be  structured  so  that  the  set  of  multiple  convolu- 

f'yn  ) 

tions  of  the  projected  correlation  function  sDj  (k)>  become  independent  of  X ^ 

and  need  to  be  calculated  only  once;  the  X ^  and  X ^  dependence  is  thus  transferred 

to  another,  simpler  part  of  the  integration.  For  correlation  functions  that  cannot 

be  transformed  to  circular  symmetry,  numerical  evaluation  is  more  costly,  since 

the  set  of  projected  correlations  |d”(1c)|  will,  in  general,  depend  on  X.  as  well  as 

on  k  and  k  separately, 
x  y 

Since  our  main  concern  is  in  studying  the  effects  of  multiple  scattering  from 
physically  realizable  turbulence,  we  will  concentrate  on  spectra  having  the  power 
law  form  given  by  Eq.  (10). 


$3d(k)  = 


*3 


where 


*3  =4" 


L.8]  1+v 


Ljl  Li/>+l) 
r(3/2)r>-i/2) 


(10) 


(id 


'N/  3d  %  3d 

is  the  normalization  constant  for  <f>  (k)  defined  by  the  condition  <P  (r=0)  =1. 

For  such  spectra,  analytic  evaluation  is  a  reasonable  possibility,  but  only 
if  V  equals  an  integer  or  half  integer.  Fcr  V ,  an  integer  or  half  integer  the 
convolution  integrals  may  be  evaluated  straightforwardly  by  contour  integration. 
However,  even  for  the  simplest  case  of  l>  =  1,  the  proliferation  of  terms  for 
higher  order  convolutions  makes  this  a  tedious  procedure.  For  the  Kolmogoroff 
spectrum  with  V  =  5/6,  the  convolution  integrations  are  best  carried  out  numer¬ 
ically.  To  accomplish  this,  we  suggested  in  I  representing  the  Kolmogoroff 
spectrum  as  a  sum  of  Gaussians.  This  allows  all  of  the  convolution  integrals  to  be 
evaluated  analytically,  yielding  a  new  sum  of  Gaussians  as  the  result.  Upon  try¬ 
ing  this  approach,  however,  we  found  that  the  series  converged  slowly,  and  that 
such  a  large  number  of  terms  was  needed  to  represent  the  Kolmogoroff  spectrum 
adequately  that  no  advantage  could  be  gained  over  direct  numerical  integration. 


To  evaluate  the  convolution  integrals  for  the  cross  section  [Eq.  (I.  54)]  most 
efficiently,  it  is  advantageous,  as  pointed  out  in  I,  to  use  the  delta  function  pro¬ 
jection  operators  appearing  in  the  projected  correlation  functions  D.  .  [Eq.  (I.  39)] 
to  integrate  out  the  z  components  of  momentum.  The  resulting  convolution  inte¬ 
grations  are  then  confined  to  the  k  ,  k  plane.  Using  th  2  results  of  the  theorem 

x  y 

proved  in  I  (Appendix  C),  all  convolutions  of  the  projected  correlation  function 
based  on  Eq.  (10)  such  as 

ii?  (k)  =  - o  D®  1  (k-k  ')  D.  <k'  )  (12) 

1  J  (27rr  1  1 


mav  be  written  in  the  k„ 


ky  plane  as 


<k^i'  k2> 


/ 


dk; 


2  ^^((kj-kj')*!^,  (kg-kJ^JDjd^^,  k^‘)  (13) 


,2  .2 


2  ns  2,2 


(27T) 


where 

=  1  +  7J2  tan2  (X.).  I?  =  L,j  /L j.  ,  (14) 

V/  -♦  _  _ 

and  where  D^(k)  is  defined  through  the  relation  [Eq.  (I.  56)  J 

D.(k)  =  D.(kx)  d  (k  •  K.)  27rKz  (15) 


as 


V  ■¥ 

Di<k  . 


)  = 


Di(kl  if. 


k2  > 


0i  (kj  kg)  = 


‘2i 


+  k^(jtiL. 


l)2+  k2 


2-jl+r 


‘-'j. 


(16) 


with 

k2.  =  47r(Lx  A.) hj.V.  (17) 

p  J  . 

The  normalization  constant  *2i  is  defined  by  the  condition  that  0  (rx=  0)  =  1 
which  yields  Eq.  (17)  when  used  with  scale  lengths  Lj  =  ,  L2  •-  L  j.  .  The 

path  dependent  attenuation  constant  a^.,  defined  through  the  relation  [Eq.  (I.  39)]  , 
is  given  by 
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(18) 


40(2  7T)liL±  cos  (*.) 


In  terms  of  dimensionless  variables  p x  =  k  x  Lj.  Eq.I13)  may  be  re-expressed  as 

%  <44  p2>  ■  I!  f  rzs  <p2-P22»£i<pl2'i2-p'22»- 

,2  7T)  (lg) 


By  rescaling  the  pj  momentum  coordinate,  we  can  rewrite  Eq.  (19)  as  a  two-dimen¬ 
sional  convolution  of  circularly  symmetric  functions.  Thus,  With  the  coordinate 
transformation 


P1  1  P1  xi*  P2  p2 


and  defining 

Kn:  Kn- 

D  (PJ - - - -  D  (p  ) 

[l+P^l  *+P  [l+p^pj 


we  have 

D”  (Pj.  ) 


(XjL 


2.n-l 

l  ) 


(20) 


(21) 


(22) 


The  multiple  convolution  products  for  the  scaled  projected  correlation  function 
will  then  be  given  by 


2 

Di  <pf> 


■  /°r‘ 


(  (Px-  P’x)  ) 


®j(p  J 


dPi 


(2  nr 


(23) 


Using  the  results  of  the  Appendix  for  convolving  two  circularly  symmetric 
functions,  Eq.  (23)  may  be  written  equivalently  as 


(25) 


•--.2  Q^MQ+P)2 


D”"1(Q2)Di(Q'2)  dO" 2 / 2 


1  7  dQ^  H  ~  /  i'  ^  ' 

(2TT)2  JQ  2  Jq*2  =  (Q  -  P)2  [(Q  +  P)2  -Q'2]  1/,2[Q2-(Q-P)2]  1/2 


or 
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=  J  QdQ  J  -  --  2  Di  <Q  )Di«?  > 


4(Q4Q*  -  P) 


(26) 


O  (2  7T) 


where  A(Q  +  Qf  -  P)  =  1  if  0,  Q',  and  P  form  the  sides  of  a  triangle  of  area  A 
and  zero  otherwise. 

By  making  the  coordinate  transformation  given  by  Eq.  (20),  integration  of 
a  general  term  such  as  Eq.  (1-60)  can  be  re-expressed  as 


dk'  E>P  (k)<<5€2(k-k'  -  K.  +  K.)  >Si  (k* ) 
JJ  (27r)6  i  ,  i  f  f 


(27) 
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-  ffff 
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The  A’,  and  dependence  of  the  integrand  has  thus  been  transferred  from 
the  arguments  of  and  ,  where  it  appeared  in  the  form  of  the  scale  factors 
SL  ^  and  Ip  to  the  term  <d€^>.  The  factors  n+^  Lj.2n  will  be  cancelled 

when  account  is  taken  of- the  terms  K  2^  *2j  appearing  in  . 

For  spectral  correlation  functions  that  do  not  lead  to  circular  symmetry  in 

the  scaled  k  -k  plane  after  integration  over  the  k  projection  delta  function, 
x  y  z 

such  as  the  factored  #Kq  of  Eq.  (1.30)  below: 
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;.'-rVNC*, 
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■thijhfcL 


BV»r"  1,(1  e 


(28) 
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[1  +k2Lx2]1  +  ,/ 

V„ 

each  Dj  (Pj»P2)  will  be  a  function  of  and  of  Pj  and  p2  separately.  While  the 
time  required  to  evaluate  a  look-up  table  of  values  for  the  non-circularly  symmet- 

V 

ric  D.  (pj,  p2)  functions  is  roughly  equal  to  that  required  for  the  symmetric  ones, 

the  storage  requirements  will  be  much  larger  since  the  circularly  symmetric 

structure  functions  D”  (P)  depend  only  on  one  variable  P  =  [  P2  +  P2  ]  *^2 

(aside  from  the  easily  factored  dependence  on  X.  stemming  from  ). 

In  the  following  sections,  we  shall  discuss  the  results  of  numerically  evaluating 

the  multiscatter  cross  section  for  4>  gi  ven  by  Eq.  (10).  We  have  omitted  the 

evaluation  of  multiscatter  cross  sections  using  the  mixed  representation  Eq.  (28) 

simply  as  a  matter  of  economy.  If  we  neglect  multiscatter  contributions,  however, 

the  attenuated  Born  cross  sections  for  the  Gaussian  [Eq.  (1. 29)],Kolomogoroff 

[Eq.  (10)], and  mixed  [Eq.  (28)]  spectral  correlation  functions  can  be  readily 

compared  and  are  shown  in  Figure  1  for  transverse  and  longitudinal  scale 

lengths  of  Lx  =  I,  Lu  =  4  and  wavelength  A.  =  1  for  A".  =  45®  .  At  the  angle  for 

specular  backscattering  X[  =  135®  ,  we  have  AKz  =  0.  Thus,  the  magnitude  of 

the  backscattering  peak  will  be  controlled  by  the  magnitude  of  AK±  •  Lx  and  the 

overall  Rayleigh  (2  7r/A)4  factor.  A  change  in  wavelength  from  A  =  1  to  A  =  1/2 

will  entail  an  approximate  25  orders  of  magnitude  decrease  in  the  Gwussian  cross 

section,  whereas  the  two  Kolomogoroff  cross  sections  are  virtually  unchanged 

due  to  the  K4  Rayleigh  prefactor  increase  balancing  the  (1+AK^  L2)  1  +  5^6 

denominator  increase.  In  the  forward  direction,  both  AK~  and  4Ki  are  zero, 

^  4 

and  the  peak  cross  section  values  are  determined  by  the  Rayleigh  K  behavior. 


4.  SLAB  GEOMETRY 

In  I,  the  symmetry  axis  of  the  irregularities  was  chosen  to  be  parallel  to  the 
face  of  the  slab.  It  is  a  straightforward  matter  to  modify  the  cross  section  calcu¬ 
lation  to  account  for  a  geometry  having  the  irregularities  aligned  at  an  arbitrary 
angle  with  respect  to  the  face  of  the  slab. 

For  the  general  case  shown  in  Figure  2,  we  start  from  the  situation  of  a 
horizontal  slab  with  irregularities  aligned  perpendicular  to  the  slab  face  and  then 
define  a  tilt  angle  a  that  defines  the  rotation  of  the  symmetry  axis  of  the  irregu¬ 
larities  in  the  xz  plane  away  from  the  z  axis,  and  a  rotation  angle  (3  that  defines 
the  rotation  of  the  projection  of  the  irregularity  in  the  xy  plane  clockwise  from  the 
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>  v 


Figure  2.  Geometry  of  Scattering  Angles  With  Respect  to  a  Fixed  Slab 
With  Rotated  Irregularities 


* 

x  axis.  The  path  length  factor  s^(x)  will  now  be  replaced  by 
z 

s'.(z)  =  -  (29) 

Icos^l 

with  z  measuring  distance  into  the  slab.  The  path -dependent  attenuation  and  the 
scattering  cross  section  can  now  be  found  by  expressing  the  wave  vector  components 


*  In  I,  the  angle  X  in  Figure  A1  was  mislabeled.  It  should  be  the  angle  ACD 

yielding  an  expression  for  the  path  length  Eq.  (1.43)  s(x)=x/  |cos(#)sin(  X )!  * 
AC.  (Thi3  will  not  entail  a  change  in  the  numerically  computed  curves  in  I 
as  4> j  was  taken  to  be  0  or  7T  for  all  cases). 


parallel  (Kg)  and  perpendicular  (Kj  2>  to  the  irregularities  in  terms  of  two  coor¬ 
dinate  rotations  through  angles  ot  and  fi . 

Thus  we  have 


Kj  =  Qj'  cos  a  +  Qg'  sin  a 

K2  =  S' 

Kg  = -Qt  "sina  +Qg'cosa  (30) 

with 


Qj'  =  Qx  cos  0  -  Qy  sin  /? 
Q2"  =  Qx  sin  &  +  Qy  cos  /? 


Qg'=Q 


z 


where  Q^.  Qzf  =  { K  sin  (0)  cos  ($),  K  sin  (0)  sin  (0),  K  cos  (0)}  are  the 

components  of  the  wave  vector  with  respect  tc  the  slab  with  rotated  irregularities, 
and  {Kj,  Kg,  Kg}  =  |k  sin  ( X )  cos  ( <f> ),  K  sin  (X)  sin  ( <t> ),  K  cos  (X)}  are  the 
components  of  or  Qf  with  respect  to  the  irregularity.  The  steps  required  to 
determine  q_q^  (<*./?)  for  the  general  case  of  scattering  from  a  horizontal  slab 

with  rotated  irregularities  are:  (1)  for  given  Q.  and  Q,,  use  Eq.  (30)  to  find  the 
corresponding  K  and  K,;  (2)  use  the  results  of  Eq.  (I.  40)  to  find  y .  [q.] 

ana  Ct  ^  LQf  J  J  (3)  according  to  whether  Q.  and  are  on  opposite  sides  or  the 
same  side  of  the  slab,  use  [along  with  5. 1  for  s'.  .(z)]Eqs.  (6)  and  (8)  or  Eqs. 

(7)  and  (9)  for  the  attenuation  and  V  factors;  (4)  determine 


Qi' 


"KtQl]Sf[<%3(aXi.  J[®1.  I1’ 


(31) 


As  an  example  of  the  effect  of  rotating  the  irregularity  axis  on  the  multiple 
scattering  cross  section,  we  take  /3  =  0  and  evaluate  °  Q^,  for  various  values 
of  Ct  for  the  Kolmogoroff  spectrum.  Assuming  that  initial  and  final  wave  vector 
are  confined  to  the  <t>  =  0  (x-z)  plane,  we  may  redefine  the  ranges  of  X^£  and 
0.  ,  to  extend  to  ±2  7T. 

If  I 

In  this  plane,  we  then  have  the  relation 
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(32) 


and 

Aff  =  a  +  d{ 

as  shown  in  Figure  3. 

For  soecular  scattering  n  -  A\  =  -  n  and  from  Eq.  (32),  this  implies 

0f  =  -2 a  -  ev  (33) 

For  6  ^  =  4?  and  for  values  of  a  =  0°  ,  30°  ,  60°  and  90°  ,  this  yields  expected 
outgoing  angles  for  spectral  scattering  of  6 ^  =  315°  ,  255“  ,  195°  ,  and  135° 
respectively. 


Figure  3.  Geometry  of  Scattering  Angles  With  Respect  to  a 
Fixed  Slab  With  Rotated  Irregularity  Axes  for  /5=  0 
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In  Figure  4,  we  have  plotted  the  total  cross  section  <rn  (a,  p  *0)  vs  0. 

for  values  of  a  =  0°  ,  30°,  60°,  and  90*  for  a  slab  with  L=  1.  rj  =  Lii/Lj.  =  4 
with  K  =  47T  for  Kolmogoroff  turbulence.  The  spectral  peaks  appear  close  to  their 
expected  values.  While  there  will  be  narrow  ranges  near  6 f  =  90°  and  270°  where 
the  cross  section  will  vanish  due  to  attenuation  in  an  effectively  infinite  slab,  the 
decrease  in  0q  ,  q  for  the  ranges  of  approximately  ±  10°  on  either  side  of 

these  grazing  angles  does  not  show  up.  This  is  caused  by  a  near  cancellation 
between  the  attenuation  of  the  Born  term  and  the  added  contribution  of  multiscatter 
terms  with  n  >  1  as  can  be  seen  in  Figure  5  where  the  two  contributions  are 
plotted  separately. 


5.  NUMERICAL  RESULTS  FOR  MULTISCATTER  CROSS  SECTIONS 

In  this  section,  we  present  numerical  results  of  multiple  scattering  cross 

sections  as  a  function  of  the  polar  angle  X  t  for  various  values  of  0  =  L, 

1  Lx 

K,  and  A-..  These  results  extend  the  calculations  of  I  to  forward  directions  as 
discussed  in  Section  2  and  to  the  case  of  the  Kolmogoroff  spectrum  as  discussed 
in  Section  3.  For  all  cases  considered,  we  assume  scattering  from  a  slab  of 
irregularities  with  the  long  axis  of  the  irregularities  parallel  to  the  slab  face 
( a  =  90°  jS  =  0 ).  In  keeping  with  the  angle  designations  of  I,  we  shall  plot  cross 
sections  vs  Af .  For  purposes  of  displaying  results,  we  define  A  f  to  be  negative 
in  the  forward  direction,  Af  =  -  |  Af  |  for  <p ^  =■  0. 

In  Figures  6  and  7,  the  logarithm  of  the  total  scattering  cross  section  a  y .,  y 
is  plotted  vs  y^  for  the  Kolmogoroff  and  Gaussian  spectral  functions 
respectively.  In  the  forward  direction,  the  Kolmogoroff  cross  section  is  seen  to 
satisfy  the  requirements  of  the  forward  scatter  approximation  better  than  the 
Gaussian  in  that  the  cross  section  falls  off  much  more  rapidly  as  |yf  -  y  .  | 
increases.  Around  the  specular  backscattering  peak,  it  is  the  Gaussian  cross 
section  that  falls  off  most  rapidly  as  a  function  of  y^.  This  behavior  can  be  readily 
understood  in  terms  of  the  derivatives  of  the  corresponding  spectral  functions. 

QX  In  4>g  =  -^-|aKx  cos  (Af)  cos  l!  -  AK„  sin  (Af)  L„2|  (34) 

n  v  =  -  K  ^K.  cos  ( X ,)  cos  Lx  -  AKh  sin  (v  )  L||2} 


[•)»] 
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SS  SECTION 
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Kolmogoroff  curve  is  approximately  four  times  that  of  the  Gaussian 


~5 


In  *K  ^2.(1  + *G 


(36) 


The  factor  in  brackets  in  the  numerators  of  Eqs.  (34)  and  (35)  can  be  rewritten  for 
the  backecattering  case  (^  =  0,  =  tt  )  as 

(37) 


[AKX  cos<Xf)L2x  ”  AKii  sin(Xf)L?|] 


=  -2K  sin 


in^j  +Xf|c°s  /xj  ~  Xf'j 


Xi  ~  Xf)  cos(Xf)L<i  -  sin 


M si 


sin(  Xf^Ltl 


For  L||  »  l/i  this  term  will  have  two  zeros,  one  near^  =  0  and  one  near  the 
specular  direction  Xf  =  Xi  • 

As  TJ  is  reduced,  the  peak  near  the  specular  direction  will  move  into  lower 

Xf  while  the  minimum  at  x  f  =  0  will  move  out.  At  some  value  of  T]( 8  <  >7  <  16 

for  the  examples  shown),  the  peak  and  minimum  will  merge  and  disappear.  This 

behavior  can  be  seen  for  both  the  Kolmogoroff  and  Gaussian  spectra.  In  the  case 

of  the  Kolmogoroff  spectrum,  the  extra  term  in  the  denominator  of  Eq.  (35) 

K  G 

implies  that  the  slope  of  In  o  will  be  smaller  than  that  for  In  a  by  2(1  +i>)/  [  1  + 
Ak\  Li  +AK||^Li|2]  which  is  approximately  2(l+5/6)/|_l  +  2K^]  o;  1/14  for  the 
parameters  of  Figure  6. 

The  equal  spacing  in  the  values  of  the  amplitudes  at  the  scattering  peaks  in 
Figures  6  and  7  reflects  the  multiplicative  factor  Lli  in  the  normalization  constants 
for  and  .  Since  the  increments  in  Ln  in  Figures  6  and  7  proceed  by  fac¬ 
tors  of  2,  the  increments  in  <r^  and  <7q  at  their  peaks  will  be  additive  in  incre  - 
meats  of  log  2. 

In  Figures  8  and  9,  we  have  plotted  the  multiple  scattering  cross  sections  for 
the  Gaussian  and  Kolmogoroff  spectra  for  increasing  values  of  K.  Fcr  the 
Gaussian  case  shown  in  Figure  8,  we  have  plotted  the  total  cross  section  for  L=l, 

T}  =  4,  and  for  K  taking  values  of  1.  5  7T,  2  7T ,  2.  5  7 t,  3  7T,  and  3.  5  7T  .  At  the 
position  of  the  forward  peak,  since  the  scattering  slab  is  relatively  thin  and 

AK  =  0,  the  K  dependence  of  the  scattering  cross  section  is  determined  by  the 

A 


Rayleigh  K  factor.  At  the  backscatter  peak,  as  pointed  out  in  1,  the  cross  section 

Thus  for 

.-23 


2  2  2  2 
is  dominated  by  the  exponential  factor  e-(AKx  Lx  +  AKu  L|i)/4. 


K  =  3.  5  it  the  value  of  the  Born  contribution  to  the  backscatter  peak  is  0. 17  x  10 

-22 

while  the  multiscatter  contribution  there  is  0. 18x10  .  Eecause  the  CFSB 

approximation  depends  on  satisfying  KL  »  1,  the  Gaussian  spectrum  offers  little 
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Figure  9.  Logarithmic  Plot  of  Total  Cross  Section  vs  Xf  for  Kolmogoroff  Tur¬ 
bulence  for  K  =  47T(+).  K  =  6 7T(x),  K  =  8lf(0)»  and  K  =  107r(  f  >  With  TJ  =  4,  L  =  10, 
6e  =  o.oi 


of  Interest  in  the  region  of  specular  scattering.  The  backscatter  cross  section 
there  for  large  KL  is  practically  zero. 

In  Figure  9.  we  have  plotted  the  total  cross  section  for  scattering  from 
Kolmogoroff  turbulence  with  V  =  4,  L  =  10,  and  for  K  values  of  4  7T «  6  it  ,  8  n. 
and  10  7T.  For  these  values  of  K  and  L,  the  forward  approximation  is  adequately 
satisfied.  In  the  forward  direction,  the  scattering  peak  is  determined  by  the 

4 

Rayleigh  K  factor  modified  by  attenuation  because  L=10  represents  a  relatively 
thick  scattering  slab.  At  the  specular  peak  near  Xf  *  40° ,  the  value  of  is  seen 
to  change  by  a  relatively  small  amount  as  a  function  of  K.  For  the  Born  contri¬ 
bution  this  behavior  is  caused  mainly  by  K^/  [l  +  A  K2  L.f  J1  +  ~  K1^2 

as  K  increases  (modified  by  attenuation  e  with  a«K  ).  The  multiscatter 

contribution  increases  roughly  as  a,  and,  for  K  =  10  n  ,  can  begin  to  dominate 
the  Born  contribution.  This  is  shown  in  Figures  10  and  11  where  we  have  plotted 
the  Born  and  multiseatter  contributions  for  K  =  10  7T  and  K  =  4  n  . 

It  should  be  pointed  out  that  the  CFSB  approximation  is  also  limited  to  domains 
where  aL  <  it  12.  Since  the  average  phase  variance  is  given  by  <.<f>  2  >  = 
aLeff  »  the  point  aL^  =  n  12  marks  the  beginning  of  the  diffusion  regime  where 
transport  theory  must  be  employed.  2  This  boundary  will  always  be  crossed  for 
scattering  at  near  grazing  angles  to  the  slab,  since  the  effective  path  length 
-*  °°  for  %  f  =  0  [0j  =  JT/2]  .  As  K  increases,  this  domain  is  bounded  by 

±  4°  ,  ±  9°  ,  ±  14°  and  ±20°  for  K  =  4jr,  6jr,  8  IT  ,  and  10  7 T  respectively 

for  the  Kolmogoroff  cross  sections  of  Figures  9,  10,  and  11. 

In  Figures  12,  13,  and  14  we  have  plotted  the  attenuated  Born  and  multiple 
scattering  contributions  for  the  Kolmogoroff  spectrum  for  V  -  4,  K  =  8  7T  ,  and 
for  L  =  0. 1,  1,  and  10.  For  the  thin  slab  case  of  Figure  12,  the  multiscatter 
contribution  is  less  than  1  percent  of  the  total  except  in  the  region  around  the 
grazing  angle  at  \  ^  ~  °°  .  Here  the  Born  term  decreases  because  of  the  attenua¬ 
tion  factor  given  by  Eq.  (6)  or  (7),  while  the  multiscatter  term  increases  due  to 
the  increased  effective  path  length.  For  the  medium  thick  slab  of  Figure  13,  this 
behavior  persists,  except  that  very  near  X[  =  °°  »  attenuation  begins  to  dominate 
the  multiscatter  term.  For  the  thick  slab  case  of  Figure  14,  the  multiscatter  term 
is  comparable  to  the  attenuated  Born  term,  and  the  domain  where  attenuation 
dominates  the  multiscatter  term  near  X[  =  0°  has  increased.  The  angular  domains 
where  diffusion  scattering  is  appropriate  are  given  by  (-2°  ,  +2  °)  for  the  case 
L  =  1,  and  (-15° ,  +15°  )  for  the  case  L  =  10. 

The  dependence  of  the  total  multiple  scattering  cross  section  on  x  f  is  shown 
in  Figure  15  for  the  case  Tj  =  4,  L  =  1,  K  =  2tt  for  incident  angles  of  X  j  =  15°  , 
45°,  75°  ,  and  85°  .  A  feature  to  be  observed  for  the  X  •  dependence,  as  pointed 
out  in  I,  is  that  the  curves  peaked  around  smaller  values  of  |  X %  I  are  broader 
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Figure  10.  Logarithmic  Plot  of  Born  (*),  Multiscattering  fy'),  and  Total  Cross 
Sections  vs  Xf  for  Kolmogoroff  Turbulence  for  K  =  107T,  T)  =  4,  L  =  10,  be  =0.01 
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Figure  15.  Logarithmic  Plot  of  Total  Cross  Sections  vs  Xf  for  Kolmogoroff  Tur¬ 
bulence  for  Xi  =  15°  (+),  X ,  =  45°  (0).  Xj  =  75°  (0).  Xt  =  85°  (4)  With  K  =  27T,  r)  =  4, 
L  =  1.0,  if  =  0.01 
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since  we  have  plotted  with  an  abscissa  proportional  tody  rather  than 

dKf  =  sint^fj)  •  K  •  d^  .  Thus,  increments  of  X  f  near  X^  =  0  yield  smaller 
2  f 

changes  in  dK^  and  therefore  broader  peaks.  The  observed  decrease  in  the 
z 

magnitude  of  the  peak  in  the  specular  backscatter  direction  where  4K||,  =0,  is 
governed  by  the  factor  1/ [l  +  A  Kx  L?.]*+l/  =  1/ [l  +  4K^  sin^  0^  Lx  ^ 
which  decreases  with  increasing  $ 

JL 


6.  SUMMARY 

We  have  presented  a  method  for  calculating  multiple  scattering  cross  sections 
from  slabs  of  arbitrary  depth  containing  a  homogeneous  distribution  of  refractive 
index  irregularities.  The  irregularities  are  assumed  to  be  aligned  along  a  sym¬ 
metry  axis  at  a  given  angle  to  the  slab  face  and  to  possess  different  correlation 
lengths  along  and  perpendicular  to  this  axis. 

Previous  calculations  that  have  been  carried  out  for  multiple  scattering  cross 
sections  have  been  limited  to  isotropic  turbulence  having  a  Gaussian  power  spect¬ 
rum  and  to  direct  backscattering.  The  extension  to  anisotropic  Kolmogoroff 
spectra  and  scattering  at  arbitrary  angles  thus  presents  a  result  that  corresponds 
to  the  general  situation  encountered  in  scattering  hf  waves  from  ionospheric 
irregularities. 

For  correlation  functions  that  depend  on  a  single  radial  argument 
.  r  2  2  o  2  2 1 

R  =  V  L  ( Xj  +  x2  )  Lx  +  Xg  L n  J  ,  projections  along  a  given  direction  in  the 

medium  can  be  changed  by  a  scale  transformation  to  circularly  symmetric  form. 

Numerical  evaluation  of  multiple  scattering  cross  sections  for  turbulent  spectra 

represented  by  such  correlation  functions  can  then  be  carried  out  relatively  simply 

and  were  evaluated  for  an  anisotropic  Kolmogoroff  spectrum. 

3  4  5 

Measurements  of  ionospheric  turbulence  spectra  *  *  tend  to  indicate  a  power 
.//  3 

law  spectrum  k  with  ft  -  2-4.  Erukhimov  et  al  further  indicate  that  the 
factored  spectral  form  of  Eq.  (28)  is  supported  by  their  data. 

We  have  noted  that  the  loss  of  circular  symmetry  in  the  plane  of  the  projected 
correlation  function  for  spectral  forms  like  Eq.  (28)  requires  more  extensive 
numerical  calculation  because  the  results  of  each  convolution  must  be  stored  as  an 
NxN  array  rather  than  as  an  Nxl  array.  If,  however,  the  initial  and  final  scatter¬ 
ing  angles  are  such  that  tan  X .  »  2(LX  /Lii  ),  the  projected  correlation 

r 

function 
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K 


v  _* 
<k) 


•  e 


-k2  tan2(A'  )Ln  2/4 


can  be  approximated  by 


Dx(k)  - 


.  -k2  tan2(A')Ln  2/4 


2.  2-1 1 +U 


[l+k^Lx 


and  the  two-dimensional  convolution  integrals  factor  into  two  one -dimensional 

convolution  integrals  of  which  the  Gaussian  part  can  be  evaluated  analytically. 

2 

Using  Erukhimov's  values  of  L  ~  0.  5  km  and  Ln  ~  10.0km  the  above 

analysis  indicates  that  for  X,  >  ~  20°  the  simpler  one-dimensional  convolu- 

x,f 

tion  calculation  should  be  adequate  and  would  constitute  an  upper  bound  to  the  full 
numerical  calculation  of  multiple  scattering  contributions. 
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Appendix  A 


To  show  the  equivalence  of  Eqs.  (20),  (21),  and  (22),  we  start  from  the 
expression  for  the  two-dimensional  convolution  of  2  arbitrary  two-dimensional 
functions  gtej.Xg),  ffXj.Xg)  whose  two-dimensional  Fourier  transforms 
g(p),  f(p)  exist. 


d(k  )  =  /  e^'*  f(x)  g(x)  dx 

-  /  dx  /  e*4*  /  e'^% 


_da_ 


_d£_ 


(2  it)2  (2  7T  )2 


(Al) 

(A2) 


1 _  f  •*■*■■*■  rv  -*■  .  -*• 

(2  7T  )2  J  6  (k-q-p)f(q)g(p)dq  dp 


->  -V  -+ 


f(q)g(k-q). 


(A3) 

(A4) 


For  functions  f(p)  and  g(p)  that  depend  only  on  the  magnitude  of  p,  this  may  be 
written  as 

•  27 r 


d(k) 


0J  <2  TT  r 


f(q2)  g(k2-2kq  cos (4>)  +  q2). 


(A  5) 


Going  back  to  Eq. 


(A2),  we  may  re-express  this  as 
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d(k) 


*  2n 

fj  fpdp  pdpqdq  fff  d<£qd0pd  <p 
eikp  cos(  ^k-  4>)  e-iqp cos  ( fpq  -  <f>  )  0-ip p  cos(<j>  .  -  <£) 


f  (q)  g(p)  /(27T)4 

oo  » 

J  pip  //•  (2  n  )3  JQ(q p)  JG(pp  )JQ(kp  )  f(q2)  g(p2 


o  (2  n  ) 


(A6) 
(A  7) 


Integrating  the  Bessel  functions  over  p ,  this  may  he  written  as 


/-X  -T-. 

dOc) 


/pdp  jqdq.  ^  s>tp2)  4(p  +  q-k) 

<2  7T  )2  A 


(A  8) 


where  4  (p  +  q  -  k)  =  1  if  p,  q,  k  form  the  sides  of  a  triangle  of  area  A  and 


k 


zero  otherwise.  Expressing  the  area  of  the  triangle  as  A  =  s(s-k)(s-p)(s-q) 

where  s  =  k+p+q.  this  may  be  re-expressed  as 
___ 


(A9) 


d(k) 


-£pdp  /qdq  f(q2)  g(p2)  4(p  +  q-k) 
(2  IT  )2  [s(s-k)(s-p)(s-q)] 


_ 4f  (p2)  g(q2) _ 

[<|H*)2-q2]1,2[q2-(p-k)2]l/2 


(A  10) 


(All) 


To  show  the  equivalence  of  this  expression  with  Eq.  (A 5)  in  another  way,  we 
let  z=q2,  a  =  (p-k)2,  b  =  (p+k)2,  and  re-express  Eq.  (All)  as 


d(k)  = 


(2  rr 


00  z=b 

"^2  y  pdp  f(p2)  2  ^  ' 


dz  g(z) 


z=a 


7  (b-z)(z-a) 


Changing  variables  to  z  =  b+z  +  Cb-a)  x  yields 

2  “2 


x=l 


d(k)  = 


l7r  )2  ^  pdp  f(p2)  / 

o 

Letting  x  =  cos  <f>  then  gives 


dx  o  o 

<2*  j  '  Jx  ^  S(P  +k  +  2pkx). 


27T 


d(k>  -  {2  7T  )2  l  pdp  f  (p2)  J  d^  g  (p2  +  k2  -  2pk  cos  <f>  ) 


= <2  7T-  )2  y*  dPx  ^2>  -  ^>2  > 


<A12) 


(A  13) 


<A14) 


which  is  Eq.(A5) 
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